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Abstract 

We present an extension of the work of D'Amato and Pastawski on electron transport in a one- 
dimensional conductor modeled by the tight binding lattice Hamiltonian and in which inelastic 
scattering is incorporated by connecting each site of the lattice to one-dimensional leads. This 
model incorporates Biittiker's original idea of dephasing probes. Here we consider finite tempera- 
tures and study both electrical and heat transport across a chain with applied chemical potential 
and temperature gradients. Our approach involves quantum Langevin equations and nonequilib- 
rium Green's functions. In the linear response limit we are able to solve the model exactly and 
obtain expressions for various transport coefficients. Standard linear response relations are shown 
to be valid. We also explicitly compute the heat dissipation and show that for wires of length 
N >> i, where £ is a coherence length scale, dissipation takes place uniformly along the wire. For 
N « £, when transport is ballistic, dissipation is mostly at the contacts. In the intermediate 
range between Ohmic and ballistic transport we find that the chemical potential profile is linear 
in the bulk with sharp jumps at the boundaries. These are explained using a simple model where 
the left and right moving electrons behave as persistent random walkers. 

PACS numbers: PACS number: 72.10.-d, 05.60.-k, 05.40.-a, 73.50.Lw 
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I. INTRODUCTION 



Inelastic scattering provides a mechanism for dissipation and decoherence in quantum 
systems. These effects are important in considering transport properties of mesoscopic 
systems. Experimental examples are numerous and include studies of transport in systems 
such as single walled carbon nanotubes^, atomic chains^'^, semiconducting heterostructures^ 
and polymer nanofibers^. In the absence of inelastic scattering, transport is either ballistic 
and we see effects such as conductance quantization^!^, or, with elastic scatterers we see 
effects of coherent scattering such as Anderson localizatiort^. In either case transport is non- 
Ohmic even when we consider very long wires. Introducing inelastic scattering necessarily 
leads to decoherence and both of the above effects (ballistic transport, localization) are 
reduced. One expects that in the limit of long wires one should get Ohmic transport^. 
Recent experiments on atomic chains^^ and Fullerene bridge have studied the effects of 
inelastic scattering and the associated local heating on quantum transport. 

The physical sources for inelastic scattering are well known and occur basically due to 
the interaction of the conducting electrons with other degrees of freedom in the system. 
For example these could arise due to electron-phonon interactions or interactions between 
conducting and non-conducting electrons^^. However the microscopic modeling of inelastic 
scattering in the context of transport is nontrivial. One of the first phenomenological models 
for dissipation was due to Biittiker— i^. In Biittiker's model one connects a point inside the 
wire to a reservoir of electrons maintained at a chemical potential fi whose value is set by 
the condition that there is no average current flow into this side reservoir. This is equivalent 
to connecting a voltage probe at some point on the wire and a nice experimental realization 
of this situation can be seen in . In Biittiker's model an electron flowing into the reservoir 
can emerge with a different phase and energy and thus one can have both decoherence and 
dissipation. 

A more detailed microscopic calculation using Biittiker's idea of incorporating inelastic 
scattering was performed by D'Amato and Pastawski^^. In their study they considered trans- 
mission across a wire modeled by the tight-binding Hamiltonian with a nearest-neighbour 
hopping parameter V. Each site on the wire is connected to electron baths which are 
themselves modeled by tight-binding Hamiltonians with hopping parameter rj. The wire is 
attached at the two ends to ideal leads with the same hopping parameter as the wire. These 
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two leads are connected to reservoirs kept at fixed chemical potentials /il and /ir for the left 
and right leads respectively. The side leads are attached to reservoirs whose chemical poten- 
tials are fixed self-consistently by imposing the condition of zero current. Using this model 
D'Amato and Pastawski analytically solved the case where the self-energy correction due to 
the side leads is pure imaginary and has the form irj and rj is small. They were able to demon- 
strate the transition from coherent to Ohmic transport. An inelastic length scale i = aV/t], 
with a as a lattice parameter, was introduced such that for wire length L << i transport was 
coherent while for L >> £ transport was Ohmic. A number of other papeiSF^^'^ have also 
shown that other models of inelastic scattering, for example due to electron-phonon (using 
side reservoirs as ensemble of harmonic oscillators to describe the heat bath) or electron- 
electron interactions, can be related to the Biittiker mechanism. Some recent papers have 
looked at electron-phonon interactions using the Keldysh nonequilibrium Green's function 
formalism combined with density-functional methods^^, tight-binding molecular-dynamics^i 
and the self-consistent Born approximation^^. An alternative mechanism for introducing 
inelastic scattering, through introduction of an imaginary potential in the Hamiltonian, has 
also been studied2^>2ii2^i2£. 

In the present paper we present an extension of the work of D'Amato and Pastawski. We 
study the case of transport of both heat and electron in the presence of inelastic scatterers 
in the form of self-consistent leads. The wire is subjected to both chemical potential and 
temperature gradients and we evaluate steady state values of both the particle and heat 
current operators. In the limit of a long wire when one is in the Ohmic regime we are 
able to obtain explicit expressions for all the linear response coefficients. It is verified that 
various linear response results such as Onsager reciprocity and the Weidemann-Franz law 
are valid. In the intermediate regime between ballistic and Ohmic transport we propose a 
simple model of right moving and left moving persistent random walkers which can explain 
much of the observed behaviour. We also perform an explicit calculation of the heat loss 
along the wire. This is a second order effect in the gradients and we show that there is 
uniform heat dissipation along the length of the wire whose value is precisely the Joule heat 
loss. For short wires we show that heat dissipation takes place primarily at the contacts. 
While heat dissipation by Biittiker probes has been discussed in 15|]27l | , we believe that this 
is the first explicitly microscopic calculation of dissipation in a quantum wire that clearly 
demonstrates Joule heat loss in the Ohmic regime and dissipation into the reservoirs in the 
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FIG. 1: A schematic description of the model. 

ballistic regime. 

The formalism used in this paper is the quantum Langevin equations approach. In two 
recent papers^^i^ it was shown how this approach can be used to derive both the Landauer 
results and more generally the nonequilibrium Green's function (NEGF) results on transport. 
Here we show how this method also works for the mulitple reservoir case and quickly leads 
to NEGF-like expressions for currents for both particle and heat. These equations are the 
starting point of our analysis. Thus apart from extending the results of [l6| we also use 
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we also consider large values of the 



a different and more general approach. Unlike 
inelasticity parameter. 

The paper is organized as follows. In sec. (ITl| we define the model and describe how the 
quantum Langevin approach can be used to get formal expressions for electron and heat 
currents in the steady state. In sec. (lllip we write the self-consistent equations and discuss 
the linear response regime. In Sec. fllVI) we solve the self-consistent equations for a long 
wire which is kept in a specified temperature gradient and evaluate the electical and heat 
current along the wire and also the heat loss into the side reservoirs. The transition from 
the ballistic to the Ohmic regime is briefly discussed in sec. (jVj). Finally we conclude with 
a discussion in sec. fIVIl). 
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II. MODEL AND GENERAL RESULTS 



We consider a one- dimensional wire modeled by the tight-binding lattice Hamiltonian. 
The wire has N sites each of which is coupled to an infinite reservoir which is itself modeled 
by a one-dimensional tight-binding system [see Fig. ([T])]. The Hamiltonian of the system 
consisting of the wire and all the reservoirs is given by 

N N 



1=1 1=1 

N-l 



where Hw = - (c^Q+i + cJ+iQ) 



1=1 

oo 



'^R = -^i E (^"^^+1 + I = 1, 2..iV 



a=l 

VU = -l'i{c^ci + clc[) l = l,2...N. (1) 

Here q and denote respectively operators on the wire and on the l^^ reservoir. The 
Hamiltonian of wire is denoted by Tiw, that of the l^^ reservoir by and the coupling 
between the wire and the l^^ site is V{y^. The coupling between the reservoirs and the wire 
is controlled by the parameters 7^. 

We briefiy indicate the steps leading to generalized quantum Langevin equations of motion 
for the wire variables. We assume that for t < to the reservoirs are disconnected from the 
wire. Each reservoir is in equilibrium at a specified temperature Ti and chemical potential 
fii. At time to we connect all the reservoirs to the wire and we are interested in the steady 
state properties of the wire. For t > to, the Heisenberg equations of motion for the wire and 
reservoirs variables are: 

Cl = ^(q_i + Q+i) + ^4 for / = 1,2...N , (2) 
cL = ^(c'.-i + cUi) for a = 2,3,...oo, / = l,2,...iV (3) 
c[ = ':^cl + '-^ci for / = l,2...iV (4) 

and we have taken cq = cn+i = 0. The equation of motion of the wire variables Eq. ([2]) 
involves the reservoir variable c^^ and we will try to eliminate this. We note that the equation 
of motion of each of the reservoirs, given by Eq. fl3lH|) . is a set of linear equations with 
an inhomogeneous part given by i'y'iCi/h . We can solve these equations of motion using 
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the single particle Green's function of the reservoirs which is given by g\t) = —i6{t)e~'^^^^/^ 
where is the single-particle Hamiltonian of the l^^ reservoir and 0{t) the Heaviside step 
function. We finally find that the solution for the boundary site on the t^^ reservoir is given 
by (for t>to) 

c\it)=^Y,9\lit-to)c'M- / dt' g[^,{t - t') ^-l c^it') (5) 

Plugging this into the equation of motion Eq. ([2]) of wire variables, we get 

ci{t) = '-^{ci.i + ci+i) - tr]i - I [ dt'Y.i{t-t')ci{t') (6) 

. I oo 

where m{t) = --j ^itit - to) c'M 

Kit) = (|)^ 

This is in the form of a generalised quantum Langevin equation where we identify rji as noise 
from the l^^ reservoir and the last term in Eq. iQ is the dissipative term. The noise depends 
on the reservoir's initial distribution which we have chosen to be an equilibrium distribution. 
The properties of the noise is written most conveniently in the frequency domain. We con- 
sider the limit to ~^ — oo. Let us define the Fourier transforms ci{uj) = (l/27r) dte^^^ciit), 
g^+ico) = JZ^dte^-'g^+it), U^) = (l/27r) J^dte^^^it) and = (yj h)' g[^,ico) . Let 

us also use the definition Ti{lj) = — /m[S^]/7r = {;~^[ / h)"^ piioj) where pi{uj) is the local density 
of states at the first site (a = 1) on the reservoir. With these definitions it is easy^S to 
show that the noise-noise correlations are given by 

{fjl{uj)f]miuj')) = Ti{uj) f{uj, Hi, Ti) 5{UJ - Uj') 5lm , (7) 

where f{u!,fi,T) = l/{exp[{!kj — ji)/kBT] + 1} is the Fermi distribution function. 

Taking Fourier transform of the equation of motion Eq. ([6]) we thus get the following 
steady state solution 

N 

ciiiv) = GU^) Vmiu;) (8) 

m=l 

where = —Z~^ 

7 

and Zim = —{uJ — S^) Sim + Sl^m-l + Sl^m+l ■ 
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As shown in 



291] G~^{uj) is basically the Green's function of the full system (wire and reser- 



voirs) and for points on the wire can be written in the form G^{uj) = [uj — Hw/h — 
where H\y is the single particle Hamiltonian of the wire while S"*", defined by its matrix 
elements S^"^ = S^^^^, is a self-energy correction arising from the interaction with the 
reservoirs. We will be interested in particle and energy currents in the system. The cor- 
responding operators are obtained by defining particle and energy density operators and 
obtaining their continuity equations^*^. The particle density is defined on sites while the 
energy density is defined on bonds. We will be interested in currents both inside the wire 
and between the wire and reservoirs. Let us define jf as the particle current between sites 
/ , / + 1 on the wire and as the energy current between the bonds (/ — 1, /) and (/, / + 1). 
Also we define as the particle current from the wire to the Z**^ reservoir and similarly 
is the energy current from the wire to the l"^^ reservoir. These are given by the following 
expectation values: 



h 



= ^{4c[-c['cn 

Using the general solution in Eq. ([8]) and the noise properties in Eq. ([7]) we can evaluate the 
above expressions and find 

N 



3i = E / Gt^G-^,^, - ) Pm {fi - U) (9) 

m=l •^^'^ 

3l = E / GUmG-^l+l - ) Pm {fl ' fm) (10) 

m=l 

= E ^ / (/^ - (11) 

m=l "^-"^ 
^ 1 

= Y.Y dcuhwTi„,{fi- U , (12) 

m=l •^-'^ 

where G^^ = G^i* and Ti^ = ^T^'^l'i^l'm PiPm\G'liJ\^ /fv" can be shown to be the transmission 
probability of a wave from the to the m*^ reservoir. 
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III. SELF-CONSISTENT DETERMINATION OF CHEMICAL POTENTIAL 
PROFILE 



We consider the case where the wire is held in a fixed temperature field specified by the 
temperatures, T/, / = 1,2...N, of the reservoirs. We will consider a small temperature 
difference and assume that the applied temperature field has the linear form 

Ti=Tl + -'—^AT , 

where AT = Tpt — T^. The chemical potentials at the ends of the wire are specified by the 
conditions fii = fJ^L and /xtv = A*/?- The N — 2 side reservoirs I = 2,3...N — 1 are included 
to simulate other degrees of freedom present in a real wire and the requirement of zero 
net particle current into these reservoirs self-consistently fixes the values of their chemical 
potentials. Thus the chemical potentials {^i} for I = 2,3...N — 1 are obtained by solving 
the following set of — 2 equations: 

f^_i = for / = 2,3,...iV-l , (13) 

with given by Eq. (fTTj) . Once the chemical potential profile is determined, we can use 
Eqs. f lQfTOj) to determine the particle and heat currents in the system while Eq. f|T2|) gives 
the heat exchange with the environment (side reservoirs). 

In general the set of equations Eq. f|T3|) are nonlinear and difficult to solve analytically. We 
will henceforth consider the low temperature and linear response regime where the applied 
chemical potential difference A/i = /x/j — ^j^l and the temperature difference AT are both 
small. More specifically we shall assume A/i << fiL,R, AT << T^ /j and ksT^^R « fJ'L,R- 
For simplicity we restrict ourselves to the following choice of parameters: 7« = 7i = 7Ar = 7 
for / = 1, 2...A^ and 7; = 7' for / = 2,3...A^ — 1. Thus all the reservoirs will have the same 
Green's function and density of states and we will use the notation g'^-^^(uj) = g^{uj) and 
pi{uj) = p{uj). 

Making Taylor expansions of the Fermi functions f{uj,fii,Ti) about the mean values /i = 
{fiL + f^R)/'^ and T = (T^ + Tr)/2, we find that, in the linear response regime, Eq. (fT3|) 
reduces to the following set of equations: 

1 TT^PT 

^--^ = E^ [^™('"'-'"™) + ^^^™(^'-^-)]=0 for / = 2,3...Ar-l ,(14) 

m=l 
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where Tim and T/^ = dTim/du are evaluated at u = ji/h. These are hnear equations in {ni} 
and are straightforward to solve numerically. We can then use Eq. ([9]) and Eq. ffTOl) to find 
the particle and heat current. The local heat current in the wire is given by = — /ijf . 
In the linear response regime we find 

N 



~ T 2^ W^lrn [fJ'l - l^m) H -F [m (-M " -imjj 



m=l 

= ^ E [ - + (^^ - ] ' (15) 

m=l 

where JF;^ = (27r277'V n^) ( G^^n'^mi+i-'^i^ im<^m/ ) P ^"^d Tim, T'lm are evaluated at = /i/ h. 
The heat loss from the wire to the reservoir can be obtained using Eq. ( |T2i) . As we shall see 
later this heat loss is a second order effect and therefore we will keep terms up to second 
order in the expansion. We then get 



3w-i 



— E 



^ '^Ira if^l ~ f^m) 7^ — %m {Tl — Tm) 



3h 

IT ^^hl _ 1 

.(16) 



+ -^m [l^l — /^m) H ~ f^m)[Ti — 1^,) H — (j] — Jm 



In the next section we will the consider the case of a long wire (A^ oo) and consider 
particle and heat transport in the presence of applied chemical potential and temperature 
gradients. Later, for an isothermal system, we will consider finite systems and discuss the 
transition from coherent to Ohmic transport. 

IV. LONG WIRE WITH APPLIED CHEMICAL POTENTIAL AND TEMPERA- 
TURE GRADIENTS 

Let us first evaluate the matrix elements Tim,{uj)- This involves p{uj) and G^^{uj). As 
discussed before p{uj) is the local density of states at the boundary site of a semi-infinite 
one-dimensional chain and is given by p{uj) = {h/{7r'j)[l — /i^a;^/(47^)]^/^ for \f}Lu\ < 27 
and zero elsewhere. For lattice points in the bulk of the wire, i.e points which are at a 
distance >> i = I/or from the boundaries of the wire we find (see Appendix [X]) = 
(-l)'+™;i e-l'-"l°/(27sinha). We now try the following solution for the self-consistent 
equations given by Eq. flTll) : 

fJ^i = f^L + j^-zj^f^ ■ (17) 



Using the fact that J2m='^oo(^ ^ ^) e"''"™-!" = 0, we see that the self-consistent equations are 
satisfied for all points I in the bulk of the wire (up to corrections which become exponentially 
small with the distance from the boundaries). Close to the boundaries the chemical potential 
variation is no longer linear. Here we focus on the limit where is very large and the linear 
solution in Eq. (ITTI) is accurate in the bulk of the wire. 

We will now use this solution to evaluate the various currents in the wire given by Eq. ( fTSjl 
and the heat loss from Eq. fll6p . We evaluate these currents at points / in the bulk of the 
wire and (since G^^ decays exponentially with distance) do not need the correct form of fii 
at the boundaries. We also find, as expected, that the currents are independent of I. They 
have the expected linear response forms: 

f = -L2iV/i-L22Vr , 

where V/x = A^/N, VT = AT/N and the various transport coeficients are given by 

1 sin^ «/ coth aji 



^11 


27vh ^ 




m= 






3 




= T L12 , 


L22 


n^klT 


3 



m] 



(18) 



7rh cosh 2aR — cos 2ai 

(20) 



^11 ! 



(21) 



where an and ai are respectively the real and imaginary parts of a and are all cal- 
culated at uj = fi/h. In deriving the above form of Lu we have used the relation 
p = 27^sinhQ;^sina//('7r7'^). In the parameter regime we are looking at, it follows that 
both Ohm's law and Fourier's law are valid, with the electrical and thermal conductivities 
given by 0" = e^Ln and k = L22- We note that Eq. fl20l) gives the Onsager reciprocity 
relation. This is usually derived within linear response theory and follows from time rever- 
sal invariance of the microscopic equations of motion. We also find from Eq. fl2Tl) that the 
Weidemann-Franz relation is satisfied. This relation states that the ratio of the thermal 
conductivity and the electrical conductivity is linearly proportional to the temperature with 
a universal constant of proportionality given by 7r^fc^/(3e^). For metals a derivation of this 
relation using semiclassical tramsport theory and within the relaxation time approximation 
can be found in {s^ . The validity of this relation requires that inelastic processes can be 
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FIG. 2: Plot of the conductivity and thermopower Qi = 3eQ / {tt'^ k'^T) (Inset) as functions of 
the Fermi level fi for different values of the inelasticity parameter 7'. They are plotted in units of 
e'^/nh and 7 respectively. 



neglected ( see discussion in [30| ). However we find that the relation continues to be valid 
in our model even though scattering is inelastic (since there is energy dissipation into the 
side reservoirs). 

From Eq f|T9|) we find that the Mott formula for the thermopower holds^i^. This is given 

by 

(3 = ^ = £Mli^. (22) 

eLii 3e crd/i 



Recently [33| have reported an interesting resonance, arising due to electron-electron inter- 
actions, observed in the thermopower as a function of the Fermi energy. We investigate if 
there are any interesting features in the dependence of Q on /i in our model. In Fig ([2]) 
we plot the conductivity and the thermopower 3eQ / {n'^ k'^T) = d{\n a) / dfi as a function of 
fi for different values of the coupling constant 7'. Surprisingly we find that for a range of 
values of the inelasticity parameter 7' there is a peak in the thermopower as a function of 
the Fermi energy. 
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Let us now look at the heat exchanges given by Eq. (fT6|) . In the long-wire limit, the 
condition of zero particle currents into the side reservoirs, Eq. f fT^ . implies that Ylim'^imil — 
m) = Ylm'^'i^ ~ ~ 0- Hence the terms linear in V/i and VT in Eq. (ITB]) vanish, and 
only the second order terms contribute significantly. Let us first consider the coefficient of 
the term containing (V/i)^ which is given by 

(23) 

Evaluating the sum we find that it is exactly equal to Lu. Determining the other terms in 
Eq. (fT6|) we find that the net heat loss per unit length (or from every bulk site) of the wire 
is given by: 

f = Lu (Vfif + iyi^)iyT) + Luiyrf . (24) 

The first term corresponds to the expected Joule heat loss in a wire and is always positive. 
The second term can be of either sign and can be identified to be the Thomson effect which 
corresponds to heat exchange that occurs in a wire (in addition to the Joule heat) when an 
electric current ffows across a temperature gradient. 

Finally we check for local thermal equilibrium in the wire. A requirement of local equi- 
librium would be that the local density ni at the point / in the nonequilibrium state should 
be the same as the density n^^'^ at the point if the entire wire was kept in equilibrium at a 
chemical potential /i; and temperature T;. It is easy to evaluate rii and and we find: 

m=l 

which in the linear response regime, gives 

,2 N 



eg 7 

m=l 



For our linear profiles of temperature and chemical potential and the form of Gim it is clear 
that, for all bulk points, the above difference vanishes (upto order 0{1/N)). Thus we see 
that the local densities are consistent with the assumption of local equilibrium. 
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FIG. 3: Plot of the chemical potential profile //j as a function of the scaled length i/N for different 
values of N and with 7' = 0.1. The points denoted by circles correspond to the approximate 
solution given in Eq. (126p . The inset shows the chemical potential profile for 7' = 1.0. 

V. CURRENT FLOW IN WIRE IN ISOTHERMAL CONDITIONS: FINITE SIZE 
EFFECTS 

In this section we look at finite length wires. We first solve Eq. f|T^ numerically to 
determine the chemical potential profile and then estimate the current in the wire using 
Eq. fllSp . We also look at the local heat dissipation at all points on the wire. 

In our numerical calculations we have chosen the parameter values /il = 1.0, = 1.1, 7 = 
1.0 and have considered different values of the dissipation strength 7' and different system 
sizes A^. In Fig. ([3]) we plot the chemical potential profile, for different system sizes, for 
a small value of the dissipation (7' = 0.1). We see that as we go to larger system sizes, 
chemical potential profile changes, from a fiat profile with large jumps at the boundaries, 
to a smooth linear profile. For a larger dissipation parameter (7' = 1.0) we see [inset of 



Fig. ([3])] that a smooth linear profi 



e is obtained even for small system sizes. The limit of 



weak dissipation was studied in 16| . Following them we find that for 7'/7 << 1 a very 
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good approximation for the transmission coefficients Tim, for any system size, is given by 

2/2/22 

%i = ' ■ (25) 

where £ = I/ur ^ 2-^"^^^. Note that, for / = 2,3, ...A^ - 1, and Ti^ are 0{i^) 
while Tim for m = 2,3, ...A^ - 1 are 0(7'''). We then find that, for £ » 1 and >> 1, 
the following chemical potential profile provides a good approximate solution of the self- 
consistent equations: 

9A 

1^1 = ^iL-S-—il-2) for / = 2,3,...iV-l (26) 
where 



2(1 + N/(.) 

Plugging in this solution into the self-consistency equations J^m'^imifJ'i — fJ'm) = with 
Tim given by Eq. ( l25l) we can explicitly verify that these are satisfied upto corrections of 
order In Fig. ([3]) we have plotted the above solution for system size = 16 and find 
an excellent agreement with the numerical result (for larger system sizes the fit agreement 
becomes better). 

The above solution leads to the following result for the current: 

= ' (27) 

where Lu is the ohmic conductivity of the wire given by Eq. f|T8l) and we have normalized 
the current such that the N ^ 00 limit gives a constant value independent of 7'. 

We have also looked at the transition from coherent to Ohmic transport for general 
values of the dissipation parameter 7'. In Fig. (jl]) we plot the scaled current jN/{A^Lii) 
as a function of system size. We find that in general for any 7'/7 < 1 the data can be 
fitted quite accurately to the form in Eq. fl27|) with i = l/an which can be interpreted as 
a coherent length scale. For 7'/7 > 1 we find that there is no coherent regime and the 
approach to the asymptotic limit has a different form. 

Persistent random walk model: It is possible to understand the various aspects of the 
intermediate regime within a simple classical Drude-like framework of right moving and left 
moving electrons moving in fixed directions but with a small probability of inter-conversion. 
We consider the case where the left reservoir is kept at a chemical potential /i + A/i and 
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the right reservoir is at /i. At the low temperatures being considered electron transport is 
basically due to the electrons close to the Fermi level and we can focus on the electrons 
within the energy gap A/i in the left reservoir. Let the density of these electrons inside the 
left reservoir be 2pL and this consists of an equal proportion of right moving electrons with 
velocity vp and left movers with velocity —vp- In the right reservoir the density of both 
left and right movers in the energy window A/i is zero. Inside the wire the presence of the 
side-reservoirs allows a right mover to be converted to a left mover with some probability. 
We now present the following random walk model to incorporate the above basic idea. The 
model consists of a lattice of N sites with a density pi of right movers and pY of left movers 
at all sites / = 1, 2...N. We impose the boundary conditions pf = Pi = Pl and Pn = P~m = 0. 
At sites I = 2, 3...A^ — 1 the particles move according to the following rules: with probability 
p a right mover at site / moves to / + 1 and, with probability 1 — p it transforms to a left 
mover and moves to site I — 1. Similarly, with probability p a left mover at site / moves to 
/ — 1 and, with probability 1 — p it transforms to a right mover and moves to / + 1. At sites 
I = 1,N, the right mover always moves to the right and the left mover moves to the left. It 
is then straightforward to write discrete time-evolution equations for the density fields pl{t) 
and pY{t). Choosing a lattice length scale a and a microscopic time scale r we obtain, in 
the continuum limit 

dp+{x,t) dp+ ^ 



dt dx 
dp'{x,t) dp 



dt dx 



+ a(p+-p-). (28) 



where v = ap/r can be identified with the Fermi velocity vp and a = {1 — p)/t gives the 
scattering rate (Note that the continuum limit requires taking a — 0, r ^ and p 1 
keeping v and a finite). We obtain a length scale v/a which we tentatively identify with the 
scattering length i introduced earlier. The boundary conditions for the above equations are 
p'^{x = 0) = Pl amd p~{x = L) = 0, where L = Na. These give the following steady state 
solution for Eq. ([2H]): 

26' 

p{x) = p+{x) + p~{x) = 2pL -6'-—x (29) 
where S' = ^^^f^ (30) 

is the density jump at the boundaries. This immediately leads to Eq. fl26l) once we note that 
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FIG. 4: Plot of the normalized current versus system size for different values of the dissipation 
constant 7'. The points denoted by circles correspond to the analytic scaling form given in Eq. (j27p . 

p{x) — 2pi (X pi — pl. The current in the wire is given by 

J = v[p^ (x) - p" (x)] = (31) 



which again leads to the result in Eq. fl27|) after we make the appropriate identifications. 

An interesting question that is often asked in the context of mesoscopic transport is: 
where is the dissipation - ? In the case of Ohmic transport, dissipation, through Joule heat 
loss, takes place in the bulk of the wire. On the other hand for coherent transport there is no 
dissipation in the bulk of the sample and the only dissipation is at the contacts (or into the 
leads). This difference between Ohmic and coherent transport can be demonstrated in our 
model by an explicit calculation of the local heat loss at all points on the wire. Using Eq. (fT6|l 
we calculate the fraction of the total heat loss that occurs at the contacts = + j^_jv 
and the bulk heat loss given by = "^f^^ Note that the total dissipation is given 

by e; =ijw-i — j^^l^ which easily follows from using the condition Yli=ijw-i — 0- The 
following table shows the contact and bulk heat losses for different system sizes and with 
7' = 0.1. In this case i ~ 200. We see clearly that for << i dissipation occurs mostly in 
the contacts to the leads while for N » i dissipation occurs in the bulk of the wire. Note 
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that the heat is eventually dissipated into the reservoirs and is possible even in a steady-state 
scenario because of the infinite size of the reservoirs. 



VI. DISCUSSION 

An interesting aspect of the present study arises if we compare it with studies of heat 
transport by phonons in oscillator chains. A big question there has been to find the necessary 
conditions on a model of interacting particles required for the validity of Fourier's law of heat 
conduct ion'^^. As a result of a large number of studies it now appears that heat conduction 
in one dimensions is anomalous and Fourier's law is not valid for momentum conserving 
models^. However there are stochastic models where one can exactly demonstrate the 
validity of Fourier's law. In one such model inelastic scattering of phonons take place by 
an exact analogue of the Biittiker probes. In this model, first proposed in [3^, and solved 
exactly recently in 37lJ38[ | , each site on a harmonic lattice is connected to a heat reservoir 



whose temperature is fixed self-consistently by the condition of zero heat current. Just as 
Fourier's law can be shown to hold in this model, here we have shown that both Fourier's 
law and Ohm's law are valid in the present tight-binding model. We have also been able to 
explicitly demonstrate local thermal equilibrium and various other linear response results. 
One other model where such a demonstration has been made in a clear way is the work by 



Larralde et al^^ on the Lorentz gas model. One other point to note is, as shown in 29[ , the 
treatments of electron and phonon transport can be done in a very similar way using the 
formalism of quantum Langevin equations and nonequilibrium Green's function. 

In this paper we have extended the calculation of D'Amato and Pastawski by study- 
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ing the finite temperature case and considering transport of botli particles and lieat in a 
tiglit-binding cliain. We liave studied botli tlie Olimic and ballistic regimes. It has been 
shown that a simple Drude-like model of persistent random walkers can explain many of the 
observed features in the intermediate regime. In the Ohmic regime we have calculated veri- 
ous thermoelectric coefficients and find that for certian values of the inelasticity parameter, 
the thermopower plotted as a function of the Fermi energy shows a peak. Finally we have 
explicitly computed heat dissipation in the wire. 

While we have only considered the linear response regime in this paper, the formalism 
described here can be used to study the nonlinear regime too. Also it can be easily used 
to study inelastic scattering effects in the tight-binding model in any dimensions and the 
reservoirs themselves can be in any dimensions. Numerical implementations to study systems 
with Anderson type of disorder and systems with externally applied magnetic fields can also 
be done readily with our approach. Finally, as pointed out in the introduction, our model 
of inelastic scattering also serves as a model for voltage probes. An important point in 
experiments involving four terminal resistance measurements on quantum wires, as in \^ for 
example, is that the voltage probe should be non-invasive. In our model the coupling to the 
probes can be tuned and thus can be used to obtain a better understanding of the role of 
probes in such experiments. Also more detailed models of the probes are easy to incorporate 
in our approach. The quantum Langevin method can be easily used for other models of the 
scattering reservoirs other than the present model where each reservoir is a one-dimensional 
wire. This would basically involve a change in the form of the self-energy correction. An 
interesting problem is an extension of the present formulation to include electron-phonon 
and electron-electron interactions. 
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APPENDIX A: EVALUATION OF GREEN'S FUNCTION 



To find the Green's function we use tlie relation = {h/'j)ZiJ^ where Z is a tridiagonal 
matrix with off-diagonal terms all equal to one. The diagonal terms are given by: 



7+ 



Zn = B{uj) = -[lo- ^g-'ioo)] for / = 2, 3...iV - 1 . 
7 ft^ 



(Al) 



The function g~^{uj) can be obtained from the green function of an isolated semi-infinite 
one-dimensional chain and, in the region of interest here < 27) is given by 



r (^) = - 

7 



27 



fc2 2\ 1/2' 



472 



(A2) 



Using standard matrix manipulations we can evaluate the inverse of Z and find 



7-1 



-1^ : for m > I 



^_^y+m£rn-lDN-l 



A 



for m < I 



(A3) 



N 



where A = AYi^i - Yi^2 

An = Det[Z] = - SAFjv-a + Yn-^ 

sinh[(/ + l)a] 



Yi 



sinh(a) 



with = - ± (— - 1)1/^ . 

2 M ' 

We will assume that the root a has been chosen such that a/? = i?e[a] > 0. Using the above 
results for the inverse of the matrix Z we find that for large the Green's function in the 
wire is given by: 



G 



Im 



27 sinh a 



-\l~m\a 



(A-e' 
{A-e- 



-(l+m-2)a 



+ e 



(A4) 
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